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Abstract 

m t 

A possible model independent test of the theoretically calculated nuclear matrix elements 
of O^/3/3-decay is proposed. The test can be accomplished if Oz^/?/?-decay of three (or more) 
nuclei is observed. The selection of the nuclei for the next generation of Ois (3 /3-decay experi- 
ments should be done taking into account considerations regarding the possibility to test the 
nuclear matrix element calculations. The test proposed allows also to check the dominance 
' of the Majorana mass mechanism of violation of the total lepton charge. 

a: 

p_; 1 Introduction 

<d : 

', The status of the problem of neutrino mixing changed drastically during the last several years: in 
.^J \ the Super-Kamiokande (SK) atmospheric neutrino |)Q, SNO solar neutrino [21 IH] and KamLAND 
reactor antineutrino (I] experiments model independent evidences of neutrino oscillations were 
^ ■ obtained. All neutrino oscillation data, except, the data of the LSND experiment |Hj 2 , can be 
described if we assume the existence of three-neutrino mixing in vacuum: 



1=1 



Here Vi{x) is the field on neutrino with mass and U is the unitary PMNS [HE] mixing matrix. 

The SK atmospheric neutrino data are best described in terms of two-neutrino — > v T oscil- 
lations. From the analysis of the data the following best-fit values of the oscillation parameters 
were found [I]: 

|Am^ 2 | = 2 ■ l(T 3 eV 2 , sin 2 2fl 23 = 1.0 (x 2 min = 170.8/170 d.o.f.)- (2) 



1 Also at: Institute of Nuclear Research and Nuclear Energy, Bulgarian Academy of Sciences, 1784 Sofia, Bulgaria. 
2 In the accelerator LSND experiment indications in favor of the transitions Pp — > P e with (Am 2 )LSND — leV 2 
were obtained. The LSND results are being tested in the MiniBooNE experiment at Fermilab 0. 
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At the 90% C.L. one has: 



1.3 • 1(T 3 < |Am 2 2 | < 3.0 • 10~ 3 eV 2 , sin 2 26 23 > 0.9. (3) 

The results of all solar neutrino experiments can be explained by u e — > v^ T transitions in 
matter. In the KamLAND experiment, v e disappearance due to transitions v e — > v^ T in vacuum 
was observed. From a global two-neutrino oscillation analysis of the solar and KamLAND data 
(performed under the assumption of CPT-invariance), the following best-fit values of the relevant 
oscillation parameters were obtained [3]: 

Am 2 21 = 7.1- 10" 5 eV 2 , tan 2 6 12 = 0.41. (4) 

In a similar 3-neutrino oscillation analysis of the solar neutrino, KamLAND and CHOOZ jU] data, 
performed in [IDj . it was found that at 90% C.L. one has: 

5.6 • 10~ 5 < Am 2 21 < 9.2 • 10" 5 eV 2 , 0.23 < sin 2 9 12 < 0.38, for sin 2 13 = 0.0, . , 

6.1 • 10^ 5 < Am 2 ! < 8.5 • 10" 5 eV 2 , 0.25 < sin 2 9 12 < 0.36, for sin 2 13 = 0.04, 

where #13 is the mixing angle limited by the reactor CHOOZ and Palo Verde experiments |9~1 lllj. 
The negative results of the CHOOZ |Hj and Palo Verde ^T] experiments are very important for 
understanding the pattern of neutrino mixing and oscillations. In these experiments no disappear- 
ance of v e was observed. From the 90% C.L. exclusion curve obtained from the analysis of the 
data of the CHOOZ experiment, the following bound can be derived 

sin 2 9 13 < 5 • 10~ 2 . (6) 

The same result was obtained in jTHj in a global 3-neutrino oscillation analysis of the solar, Kam- 
LAND and CHOOZ data with |Am| 2 | taken to lie in the interval eq. (J3J). 

There are two general theoretical possibilities for the fields of neutrinos with definite masses 
Vi(x) (see, e.g., [12]): 

1. If the total lepton charge L = L e + L M + L T is conserved, Pi(x) are Dirac fields of neutrinos 
Vi (L=l) and antineutrinos i>i (L=-l). 

2. If there are no conserved lepton charges, Vi{x) are Majorana fields which satisfy the condition 

^(x) = Cuf(x) = u t (x), (7) 
C being the charge conjugation matrix, and neutrinos i/j are Majorana particles. 

The solution of the problem of the nature of the neutrinos with definite mass - Dirac or 
Majorana, will be of fundamental importance for the understanding of the origin of small neutrino 
masses and of the pattern of neutrino mixing. 

The investigation of neutrino oscillations does not permit to determine the nature of massive 
neutrinos [T3*| IT4] . Processes in which the total lepton charge L is not conserved and changes by 
two units, must be studied for that purpose. These processes are allowed if the massive neutrinos 
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are Majorana particles. Experiments searching for neutrinoless double (3- decay of even-even nuclei 

(see, e.g., Hang nana usi eh cni), 

(A,Z) -> (A,Z + 2) +e~ +e". (8) 

have the highest sensitivity to the nonconservation of the total lepton charge and to Majorana 
neutrino masses. The matrix element of 0v/3(3 -decay is proportional to the effective Majorana 
mass (see, e.g., [12122]): 

\m ee \ = \ mi \U cl \ 2 + m 2 \U e2 \ 2 e ia21 + m 3 \U e3 \ 2 e ia31 \ , (9) 

where U e j, j = 1,2,3, are the elements of the first row of the PMNS neutrino mixing matrix U, 
rrij > is the mass of the Majorana neutrino Uj, and a 2 \ and 0131 are two Majorana CP- violating 
phases [T31I2B]- One can express [21] (see also, e.g., [221122]) the two heavier neutrino masses and 
the elements \U e j \ in \m ee \ in terms of the lightest neutrino mass, Am^, Am^, and of 612 and #13 , 
respectively. 

The results of a large number of experiments searching for Oz//5/5 -decay are available at present 
(see, e.g., [2"U1 l2"o1 |2T| I28|). In the Heidelberg-Moscow experiment j2H] the most stringent lower 
bound on the half-life of vfi /3-decay of 76 Ge has been obtained 3 : 

Ty 2 > 1.9 • 10 25 y (90% C.L.). (10) 

In the cryogenic experiment CUORICINO the following lower bound on the v(3 /3-decay half-life 
of 130 Te was recently reported |3~3j : 

T^^-lO^y (11) 

Taking into account the result of different calculations of the relevant nuclear matrix elements, 
the following upper bounds on the effective Majorana mass \m ee \ can be inferred from the limits 
(PI) and HU): 

|m ee | < (0.3 - 1.2) eV, |m ee | < (0.3 - 1.7) eV. (12) 

The NEM03 experiment searching for vf3 /3-decay of a number of different nuclei ( 100 Mo, 82 Se, 
etc.) and aiming at a precision of |m ee | ~ 0.1 eV, is successfully taking data at present 

Many new projects of experiments searching for neutrinoless double /3-decay of 76 Ge, 136 Xe, 
130 Te, 100 Mo and other nuclei, are under research and development at present (see, e.g., |19| l2"Ul 
123 ES] )• The goal of the future experiments is to reach a sensitivity 

\m ee \ ~ few x 10~ 2 eV. (13) 

In order to obtain information about the effective Majorana mass \m ee \ from the results of the 
0^/?/9-decay experiments, the corresponding 0i//?/9-decay nuclear matrix elements must be known. 
A large number of calculations of the nuclear matrix elements of the 0i//?/?-decay exist in the 
literature (see, e.g., [TBI HZ1 HH] ) • The results of different calculations differ by a factor three or 
more. We will propose here a possible test of the nuclear matrix elements calculations. It can be 
accomplished if O^/3/3-decay of several nuclei will be observed in the future experiments. 

indications for O^/3/3-decay of 76 Ge with a rate corresponding to 0.11 eV < |m ee | < 0.56 eV (95% C.L.), are 
claimed to have been obtained in |3Uj. The results announced in j3U] have been criticized in [3T]. Even stronger 
evidence has been reported recently in |32|. These claims will be checked in the currently running and future 
v(3 /3-decay experiments. However, it may take a very long time before comprehensive checks could be completed. 
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2 The effective Majorana mass 



The neutrino oscillation data allow to predict the possible ranges of values of the effective Majorana 
mass. The prediction depend strongly on the character of the neutrino mass spectrum and on the 
value of the lightest neutrino mass (see, e.g., [22 ESI EEl EZ1 EHj)- We will summarize here briefly 
the main results for the three possible types of neutrino mass spectrum jSHl HO] . 

1. Normal hierarchical neutrino mass spectrum: 

nil <^ m 2 ^ m 3 (14) 

For the effective Majorana mass we have in this case the following upper bound 



\m ee \ < (^sin 2 12 y Amlj + sin 2 9 13 ^ Am 2 2 J . (15) 

Using the 90% C.L. ranges (j2J) and (j5J) of the oscillation parameters and the CHOOZ bound 
for the effective Majorana mass one finds 



|m ee ] < 5.5 • 10~ 3 eV, (16) 

This bound is significantly smaller than the expected sensitivity of the future 0i//9/3-decay 
experiments. The observation of the 0^/?/3-decay in the next generation of experiments might 
exclude normal hierarchical neutrino mass spectrum. 

2. Inverted hierarchical neutrino mass spectrum: 

m 3 <C nil < m 2- (17) 

For the effective Majorana mass we have in this case 



|m ee | ~ J |Am| 2 | (1 - sin 2 2^12 sin 2 a) 2, (18) 
where a = a 2 i is the Majorana CP— violating phase. From eq. (|18|) we obtain the range 



Amfal cos 20i2 < \m ee \ < y|Am| 2 |, (19) 

where the upper and lower bounds corresponds to the case of CP conservation and the same 
and opposite CP-parities of neutrinos vi and z/ 2 . Using the 90% C.L. allowed values of the 
parameters, eqs. © and ©, for the effective Majorana mass one finds [39J: 

10- 2 eV<|m cc |<5.5-Hr 2 eV (20) 

Thus, if the neutrino mass spectrum is of the inverted hierarchical type and the massive 
neutrino are Majorana particles, Oz//3/3-decay can be observed in the experiments of next 
generation. 
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3. Quasi- degenerate neutrino mass spectrum: 

mi ~ m 2 — m 3 , m\ 2 3 3> | Amg 2 | . (21) 

The effective Majorana mass is given in this case by eq. (j!8)l in which a/| Amjy is replaced 
by m min , where m min is the lightest neutrino mass. For the effective Majorana mass we have 
the range [3~§] 

0.22 m min < \m ee \ < m min . (22) 

In the case of quasi-degenerate neutrino mass spectrum, the effective Majorana mass depends 
essentially on two parameters: the lightest neutrino mass m m j n and the CP-violating param- 
eter sin 2 a. From the measurement of |m ee |, the following range for the lightest neutrino 
mass can be obtained: 

\m ee \ < m min < 4.6 \m ee \ . (23) 

If the lightest neutrino mass m m i n will be determined in the tritium (3— decay experiment 
KATRIN [UJ which is under preparation at present, or from cosmological and astrophysical 
observations (see, e.g., |42j), the data of 0^/3 j3 experiments can be used to get information 
about the Majorana CP- violating phases (313 E2J IH3 E! • 

Neutrinoless double /3-decay is a unique process. The observation of this process would be a 
proof that the total lepton charge is not conserved and massive neutrinos vi are Majorana particles. 
As we have seen in this Section, the precise measurement of the parameter \m ee \ would allow to 
draw important conclusions about the character of neutrino mass spectrum, the lightest neutrino 
mass and the CP-violation associated with the Majorana neutrinos (for a more detailed discussion 
see, e.g., [13]). However, from the data of the Oi//3/3-decay experiments only the product of \m ee \ 
and the corresponding nuclear matrix element can be determined. In the next section we will 
briefly discuss the problem of the calculation of the nuclear matrix elements of Ozy/3/5 — decay. 

3 The Problem of Calculation of Nuclear Matrix Elements 

If neutrinoless double /3-decay is due to the Majorana neutrino mixing Q only, it proceeds via 
exchange of a virtual neutrino and is a process of second order in the Fermi constant Gf- The 
nuclear matrix element (NME) of the 0z//3/3-decay of a given even-even nuclei cannot be related to 
other observables and has to be calculated. The calculation of NME is a complicated problem. One 
of the problems of the calculations is connected with a large number of states of the intermediate 
odd-odd nuclei, which are important due to relatively large average momentum of the virtual 
neutrino. 

Many calculations of NME exist in literature (see, e.g., the review articles [Tol IT71 IT^l 
Two basic methods of calculations of NME are used at present: nuclear shell model (NSM) and 
quasiparticle random phase approximation (QRPA). 

The nuclear shell model is attractive from physical point of view: there are many spectroscopic 
data in favor of shell structure of nuclei (spins and parities of nuclei, binding energies of magic 
nuclei, etc.) |47j . However, only rather limited set of one-particles states of valent nucleons can 
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Table 1: Half-life of the 0^/3/3-decay for \m ee \ = 5-10 2 eV. The nuclear matrix elements were 
taken from the compilation in ref. |46j . 



Nucleus 




1 l/2 


years 




76 Ge 


1.4 


10 27 


- 1.5 


10 29 


100 Mo 


1.7 


10 26 


-5.9 


1Q 30 


130 Te 


7.7 


10 26 


-3.4 


10 27 


136 Xe 


2.7 


10 27 


- 1.7 


10 28 



be taken into account because of practical computational reasons. It is difficult to estimate the 
accuracy of the shell model calculations. 

The most popular method of calculation of the NME of 0uf3(3 -decay is QRPA ^BJ This 
method allows to use as a basis a large number of one-particle states and to take into account all 
intermediate states. Important parameters of QRPA are the constant of particle-hole interaction, 
g p h, and the constant of particle-particle interaction, g pp . The constant g p h can be fixed from a fit 
of the energy of the giant Gamov- Teller resonance. The constant g pp is a free parameter. 

There are many models based on the QRPA approach (see, reviews |lbllT7|ITm i46j). The results 
of the calculations of NME of 0z//3/3 -decay performed by different authors differ quite significantly. 
The variety of results of the calculations is illustrated by Table 1 in which ranges of the values of 
the half-life of the O^/3/3-decay of different nuclei are presented for \m ee \ = 5 • 10 _2 eV. 

Recently, in the framework of QRPA, a new procedure of calculation of the NME was proposed 
|48j . For a fixed value of the constant g p h and the values of the parameter g pp , determined from 
the measured 2z//3/3-decay half-life, the 0z//3/3-decay nuclear matrix elements of several nuclei were 
calculated. In [IB] results were derived for three different sets of the one-particle states and 
for three different nucleon-nucleon potentials, and it was shown that the NME of 0i//3/3-decay 
depends weakly on the number of one-particle states and on the nucleon-nucleon potential used 
in the calculations (the NME for each of the nuclei thus calculated varies by not more than 10%). 

Another approach was proposed in In this paper all parameters of the QRPA model, 

including the g pp constant, were fixed from data on /3-decay of nuclei which are close to the 
even-even nuclei of interest for the Ou/3/3 -decay study. 

In spite of the recent progress in the calculation of NME, it is not possible at present to estimate 
the real accuracy of the calculations. It is important to find a possibility to check the calculations 
of NME by a direct comparison with experimental data jUJ]- Such a possibility will be discussed 
in the next Section. 

4 Possible Test of the NME Calculations 

If the Majorana neutrino mixing (|1J) is the mechanism of Oz//5/5-decay, the matrix element of the 
process has the following general form (see, e.g., [TB^ IT2~]) 

(f\(S - l)|z) = Nm ee M°»(A, Z) 5(E f - E t ). (24) 
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Here N is a product of known factors and M 0u (A, Z) is the nuclear matrix element of interest. 
The neutrino masses enter into the matrix element (f\{S — through the effective Majorana 
mass m ee given by and the neutrino propagator (q 2 — m 2 ) -1 which is included in M°"(A, Z), q 
being the momentum of the virtual neutrino. For small neutrino masses (smaller than the binding 
energy of nucleons in nuclei ~ 10 MeV), the neutrino masses in the propagator can be safely 
neglected. The matrix element M 0v {A 1 Z) depends in this case only on the nuclear properties and 
strong interaction. 

The half- life of the Qv/3/3 is given by the expression 4 : 

' \m ee \ 2 \M »(A,Z)\ 2 G°»(E ,Z), (25) 



T^ 2 (A,Z) 

where G 0u (E , Z) is known phase-space factor (E is the energy release). If we use a model M of 
the calculation of NME we have 

\m ee \ 2 M (A, Z) = z) \M°f(A, Z)\ 2 G°»(E , Z) " (26) 

Let us assume that neutrinoless double /3-decay of several nuclei is observed. The effective 
Majorana mass \m ee \ cannot depend on the parent nucleus. Thus, if the light Majorana neutrino 
exchange is the dominant mechanism of 0z//3/3-decay, the model M of the calculation of the nuclear 
matrix elements can be correct only if the relations 

\m ee \ 2 M (A u Zi) ~ \m ee \ 2 M (A 2 , Z 2 ) = ... (27) 

hold, where \m ee \j^(Aj,Zj) is the value of \m ee \ 2 obtained from the 0^/3/3-decay half-life of the 
nucleus (A,-, Zj) using the model M. 

Consider different models of calculation of NME. From eq. (|2~5j) it follows that for a given 
parent nucleus the product \m ee \ 2 M (A, Z) \ \M$J(A, Z)\ 2 does not depend on the model. Thus, for 
different models and the same nucleus we have 

\m ee \ 2 Mi (A, Z)\ \M^(A,Z)\ 2 = \m ee \ 2 M2 (A, Z)\ \M°£(A, Z)\ 2 = ... (28) 

For two different models we have the relation 



\m ee \ 2 U2 (A,Z) = V M ^(A,Z) \m ee \ 2 M2 (A,Z) (29) 



where 



\M° M \{A,Z)t 

71 {AZ) -\Mir 2 (Az)r m 



4 It follows from this expression that the relative accuracy of determination of the parameter \m ee \ (for any value 
of NME) is two times better that the relative accuracy of the measurement of the half-life of O^/3/3-decay: 

AKe] _ 1 AT l/2 

\m ee \ 2 T°fi 
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Table 2: The parameter r) Mi ' Mk (A, Z), determined by eq. ([30jh for the nuclear matrix elements of 
Qup/3- decay, calculated in ref. [ED] (Mi), in ref. gHj (M 2 ) and in ref. [H] (M 3 ). 



Nucleus 








vy Ge 


0.37 


0.19 


1.93 


82 Se 




0.38 




100 Mo 






6.56 


130 Te 


0.74 


0.10 


7.32 


136 Xe 


0.53 


0.02 


22.42 



In Table 2 we present the values of the coefficient rj(A, Z) for the case of the matrix elements 
calculated in [50] (NSM), in gH] (QRPA ), and in 00] (QRPA, different model). 

We see from Table 2 that the coefficient rj(A, Z) depends rather strongly on (A, Z). This means 
that if for one model the relation (}2~Tj) is satisfied, other models, in principle, can be excluded. 
However, the observation of O^/3/3-decay of only two nuclei might not allow to distinguish between 
different models. For example, the observation of the 0z//9/?-decay of 100 Mo and 130 Te will not allow 
to distinguish the QRPA models of ref. [1H] and of ref. jlHl (the difference between the values of 
the coefficient rj(A, Z) for these two nuclei is about 10%). The values of the effective Majorana 
mass which can be obtained from the observation of the 0z//3/3-decay of these two nuclei, if one 
uses the models of ref. [48] and of ref. [IE], will differ by a factor of ~ 2.6. The observation of 
neutrinoless double /?-decay of at least three nuclei would be an important tool in the solution of 
the problem of NME. Table 2 suggests that the observation of the Qv/3/3- decay of 76 Ge, 130 Te and 
136 Xe would solve the problem.. 

If relations (|27J) are satisfied for some model M, this would mean that the corresponding value 
of the effective Majorana mass \m ee \M can be different from the true value |m ee |o by a constant 
factor: 

\m ee \ M = (3\m ee \ . (31) 

From this relation it follows that 

\M° A ?(A,Z)\ M = ±\Ml?(A,Z)\ , (32) 

where \M$f(A, Z)\ Q is the true value of the NME (which, of course, we do not know). For the test 
we are proposing it is important that nuclear matrix elements depend rather strongly on (A, Z) . 

It looks quite improbable that relation f!32l) with one and the same constant (3 would be valid 
for nuclei with different properties and different NME, and especially for three different nuclei. 
Thus, if relation (|2*7jl is satisfied for some model of calculations of nuclear matrix elements, the 
corresponding value of the effective Majorana mass would most likely be the true value. 

One last remark. We have assumed that the mechanism of Oz//3/3-decay is the Majorana neutrino 
mixing, eq. There exist, however, many other mechanisms of nonconservation of the total 

lepton charge and Oi//3/3-decay, like SUSY with violation of R— parity, etc. (see, e.g., the review 
articles These additional mechanisms modify the neutron-proton-electron vertexes and 
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include exchange not only of light neutrino but also of heavy particles (see, e.g., |KT|). Because 
the nuclear matrix elements of different operators have different A, Z dependence, in the case of 
additional mechanisms the factorization property of 0uj3j3— decay matrix elements, eq. ([24)1. will 
not be valid and relation ([27)1 . in general, will not be satisfied. Thus, if relation (|2T|) will be 
found to hold for a given model, this also will be a strong indication in favor of dominance of the 
Majorana neutrino mixing mechanism of lepton charge non-conservation. 

5 Conclusions 

After the discovery of neutrino masses and neutrino mixing the problem of the nature of neutrinos 
with definite masses Ui became one of the fundamental problems in the studies of neutrino mixing. 
The determination of the nature of massive neutrinos z/j will have a profound implications for the 
understanding of the mechanism of generation of neutrino masses and mixing. The measurement 
of the effective Majorana mass would allow to obtain an important and unique information about 
the character of the neutrino mass spectrum, the lightest neutrino mass and possibly on the 
Majorana CP-violating phases. 

From the measured half-life of Oi//3/3-decay only the product of the effective Majorana mass and 
the corresponding nuclear matrix element can be obtained. The results of the different existing 
calculations of the Oi//3/3-decay nuclear matrix elements vary significantly (see Table 1). The im- 
provement of the calculations of the nuclear matrix elements is a very important and challenging 
problem. We have discussed here a possible method which could allow to test the models of calcu- 
lation of NME of the 0^/9/9-decay via comparison of the results of calculations with experimental 
data. The method is based on the factorization property of the matrix element of 0^/?/?-decay and 
requires the observation of the 0i//9/9-decay of several nuclei. 

The nuclear matrix elements of the 0^/9/9-decay cannot be related to other observables. We do 
not see at present any alternative possibility to confront the results of the NME calculations with 
experimental data in a model independent way. 

New experiments searching for Oz//5/3-decay of 130 Te, 76 Ge, 136 Xe, 100 Mo and other nuclei are in 
preparation at present. From the point of view of the problem of NME, it is very desirable that 
these projects will be realized for at least three different nuclei. The selection of the nuclei should 
be done taking into account also the considerations discussed above regarding the possibility to test 
the nuclear matrix element calculations. 
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